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ABSTRACT 

We present an analysis of the Vela X-1 "off-states" based on Siizakii observations taken in June 2008. Defined as states in which the 
flux suddenly decreases below the instrumental sensitivity, these "off-states" have been inteipreted by several authors as the onset of 
the "propeller regime". For the first time ever, however, we find that the source does not turn off and, although the flux drops by a 
factor of 20 during the three recorded "off-states", pulsations are still observed. The spectrum and the pulse profiles of the "off-states" 
are also presented. We discuss our findings in the framework of the "gated accretion" scenario and conclude that most likely the 
residual flux is due to the accretion of matter leaking through the magnetosphere by means of Kelvin-Helmholz instabilities (KHI). 
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1. Introduction 

Vela X-1 is a persistently active high-mass X-ray binary sys- 
tem (HMXB) consis ting of a massive neutron star (I.88M0, 
lOuaintrell et"al1l2003l) and a BO.SIb type su per giant HP 77581 
with mass of ~ 23Mo and radius of ~ 30/?© dvan Kerkwiik et all 
Il995h . The average X-ray luminosity of the pulsar, Lx ~ 4 x 
10^^ ergs s"', is explained well by the mass-loss rate of the opti- 
cal companion of ~ l O'^Mpyr'' (iNagase et al.lll986h . assuming 
a simple wind model dCastor et al.lll975|) and the observed termi- 
nal wind velocity of ~ llOOkms ' dWatanabe et alj|2006[) . The 
neutron star is ecli psed by the optical comp anion every orbital 
cycle of ~ 8 .964 d (Ivan Kerkwiik et ai][ 1995b. The spin period of 
the neutron star, Ps ~ 283.5 s ( Rappaportll 1 975h . has remained 
almost constant since this discovery. The X-ray spectrum of 
Vela X-1 is described well by a cutoff power law dNagase et al.l 
Il986l) . Cyclotr on resonance scatteriiig fea tures have been re- 
orted at ~ 25 dMakishima & Miharalll992h and - 50-55 keV 
Kendziorra et al.l fl 9921) . although the feature at 25keV is still 



debated (IOrlandinill2006l) . 



The source is strongly variable with an average X-ray lumi- 
nosity of ~ 4 x lO^'^ergs"' (assuming a distance of ~ 2kpc, 
lNagasell989l) . Aside from the usual flaring activity similar to the 
one observed in other wind-accreting pulsars, abrupt "off-states", 
in which the source becomes undetectable for several pulse peri- 
ods, and "giant flares" characterized by an increase in the flux up 
to a factor of 20 have been observed (iKrevkenbohm et ani2008l) . 
Both phenomena are not unusual fo r Vela X-1, although prob - 
ably neither has a periodic nature (IKrevkenbohm et aP l2008l) . 
The giant flares of Vela X-1 have been compared to those seen 
in Su per Fast X-ray Transients (SFXTs, Walter & Zurita Heras] 
I2007h . and may be connected to the clumpy structure of the wind. 
In this scenario the wind velocity and the density fluctuations 
are responsible for the luminosity swings: flares would be asso- 
ciated with the accretion of particularly massive wind clumps, 
whereas "off-st ates" can be explained with the onset of the cen- 
trifugal barrier 

dm arionov & Sunvaev 197 53, trigger ed by a drop 
in wind density /velocity ( Krevkenbohm et al.ll2008l) . 



An alternative "gated accretion" scena rio to explain the flar- 
ing ac tivit y of SFXTs was p roposed by iGrebenev & SunvaevI 
d2007h and lBozzo et all d2008l) . who suggested that the observed 
luminosity swings may be associated to the transition between 
different accretion regimes, i.e. to the different ways the plasma 
enters the magnetosphere. It is the interaction of the rotating neu- 
tron star's magnetosphere with the plasma that mediates the ac- 
cretion rate. When the magnetospheric boundary is strong, the 
matter accumulates close to the magnetospheric boundary. If the 
growth of the magnetospheric instabilities renders the barrier 
transparent, the accumulated matter can accrete, thereby produc- 
ing a brigh t flare. Observation al ev idence for this sc enario was 
reported bv lBozzo et al ] (120081) and lGrebenevl(l2010h . 

In this work we report on the analysis of the June 2008 
Suzaku observation of Vela X-1. We focus on the three "off- 
states" detected in the observation. To our knowledge it is the 
first time that this type of activity has been observed with an in- 
strument sensitive enough to constrain the flux in the "off-state", 
to detect pulsations, and to perform the spectral and pulse profile 
analysis. 

We conclude that the accretion still proceeds during the "off- 
states", although at a much lower rate. We discus s our result in 
the framework of the "g ated accretion" scenario (iBurnard et al.l 
119831: iBozzo et al I I2OO8I) mentioned above. The abrupt decrease 
in the accretion rate during the "off-states" may then be ex- 
plained if the magnetospheric boundary becomes stable with 
respect to Rayleigh-Taylor instability, but some matter still 
leaks because the magnetosphere is still unstable with respect 
to Kelvin-Helmholz instabilities (KHI). 



2. Data analysis and results 

The observation we rely on is a ~ lOOks long Suzaku obser- 
vation (ID 403045010), performed on June 17-18, 2008, about 
1.6 d after the eclipse and close to the periastron passage of the 
source (orbital phase ~ 0-0.16). The data was reduced using 
the HEADAS 6.9 with CALDB version 20100812. 
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Timing analysis. To improve the quality of the statistics, data 
from all three XIS units were combined. The lightcurve of the 
observation in the range 0.4-12keV is presented in Fig.[T] Three 



served. As reported in literature ( 


Inoue et al.ll984llLapshov et all 


119921 iKrevkenbohmet alj 11999 


, 20081), during the "off-states" 



the sources' flux drops abruptly and recovers after several pulse 
periods. Marginal evidence o f a residual pulsed emission was re- 
ported bv llnoue et al.l ( 1 19841) based on Tenma data. With the un- 
precedented sensitivity of Suzaku not only we do unambiguously 
confirm these findings but we can also study the "off-states" in 
detail. 



Using the phase-connection technique ("Stauber t et al. 2009h 



we 



and assuming the ephemeris by Kreykenbohm et al 
determined the pulse period to be P, = 283.473(4). All uncer- 
tainties quoted are at lo" confidence level unless stated otherwise. 
The marginal evidence of spin-up is not statistically significant 
and may be attributed to the uncertainty of the orbital parameters. 
No change in the pulse period in "off-states" could be measured. 
Based on the obtained timing solution, we constructed energy- 
resolved pulse profiles for the entire observation and for the "off- 
states" by combining all three episodes. As can be seen in Fig.|2l 
the pulse profiles vary significantly with both energy and lumi- 
nosity. Very remarkable is the change around phase ~ 0.75 at 
hard energies (20-60 keV) between the normal and the "off-state" 
profile. 

We also investigated the flux dist r ibutio n of the source fol- 
lowing the approach of iFiirst et aP (l2010l) . who shows that 
the flux distribution is approximately lognormal based on 
INTEGRAL data. Although the lognormal distribution generally 
describes our data, an excess appears at low count rates (Fig.[3]l. 
This excess is due to the "off-states". The flux distiibution of the 
"off-states" is still approximately lognormal (Fig. |3]l but differs 
considerably from the distribution of the rest of the lightcurve. It 
is the "off-states" component that mostly contributes to the low- 
countrate flank seen in the overall flux histogram. Implications 
of this finding are discussed below. 



Spectral analysis We first analyzed the average spectrum of the 
entire observation to establish a baseline for the analysis of the 
"off-state" data. Several phenomenological continua, based on 
the models reported in literature for Vela X-1, were used to fit 
the average spectrum of the source. None of these models was 
able to desciibe the broad-band 0.4-70 ke V spectrum. In partic- 
ular, the spectrum below 5 keV is poorly described by cut-off 
power-law models. 

We used two components to model the continuum, c ombin- 
ing a Comptonization model by iTitarchuk et al.l (Il996h and a 
power law. Photoelectric absorption at lower energies, a num- 
ber of emission lin es, and an iron absorption edge at ~ 7.26 
(iNagase et al.l 1 19861) were also necessary to fit our data. Two 
CRSF harmonics were also required by the fit and were modeled 
using a multiplicative Gaussian profile. The best-fit parameters 
are summarized in Table [T] The quality of the statistics is sig- 
nificantly lower for the "off-state" data, so a simpler model was 
used in this case: an absorbed Comptonization model with the 
addition of an iron absorption edge. No other spectral features 
were required by the fit. The results of the best fit are presented in 
Table [Tj and the best-fit spectra and the fit residuals are shown in 
Fig.E] The average absorption-corrected flux in the 0.4-70 keV 
energy range was ~ 3.8 x 10"^ erg s"' for the complete observa- 
tion and ~ 5 X 10"'" erg s"' for the "off-state" spectrum. Results 



of a more detailed spectral analysis, including phase-resolved 
spectra, will be published elsewhere. 

It is interesting to note that the "off-state" spectrum differs 
considerably from the spectrum observ ed during the ecUpses 
when it is dominated by emission-lines (IWatanabe et al.ll2006l) 
originating in the surrounding plasma illuminated by the X-rays 
emitted by the eclipsed pulsar 

3. Interpretation and discussion 

First, we would like to summarize the observed properties of the 
source in its "off-state" 

- The flux drops by a factor of ten or more on a timescale com- 
parable to the pulse period. The source remains "off" for sev- 
eral pulse periods, and then the flux is restored to the previ- 
ous level on the same short timescale. 

- Vela X-1 is observed to pulsate in the "off-state". No pulse 
frequency change has been detected in the "off-state". 

- A drastic change in the shape of the pulse profile shape ap- 
pears at high energies. The narrow dip at pulse phase ~ 0.75, 
seen in the profiles of the normal state at all energy ranges 
as well as in the "off-state" profiles at lower energies, is sub- 
stituted by a prominent peak in the "off-state" profile at hard 
energies (20-60) keV. This leads to a significant increase in 
the fraction of the pulsed emission in the hard energy range 
for the "off-states". 

- The flux distributions of the normal and "off-states" are sig- 
nificantly different. The overall distribution is composed of 
two approximately lognormal peaks. 

- Although poorly constrained at high energies, the spectrum 
changes significantly during the "off-state". The temperature 
of the Comptonizing medium increases, whereas the optical 
depth decreases. No CRSF is required by the data, although 
the statistics at high energies does not allow us to rule out 
the presence of a CRSF completely, especially if the CRSF 
is shifted to higher energies. The absorption corrected flux in 
0.4-70keV energy range is ~ 5 x 10""^ergscm"^ s"', which 
coiTesponds to a luminosity of ~ 2.4 x lO^'^ergss"' for a 
distance of 2 kpc. 

The observed pulsations, the luminosity and the hard spectrum 
of the "off-states" can only be explained if the emission is pow- 
ered by the accretion of plasma onto the magnetized neutron 
star. The absence of emission lines in the "off-state" spectrum 
strongly suggests that the source is not eclipsed, but instead ex- 
hibits an intrinsic drop in luminosity, hence in the accretion rate. 
Th e timescale of the state tran sition makes it difficult, as argued 
by iKrevkenbohm et aP (l2008h . to explain the onset of the "off- 
states" with a sudden decrease in wind density and/or velocity, 
and suggests a magnetospheric origin of the state transition. 

This agrees with the observed flux distribution. If the "off- 
states" were due to drops in the wind density, one would expect 
them to contribute to the lower-flux tail of the normal flux dis- 
tribution. On the contrary, they form a distinct low-flux peak as 
observed in Fig.|3] The lognormal flux dis tribution is most likely 
caused, as discussed bv lFiirst et al.l (1201 Ol) . by the "grinding" of 
a clumpy wind by the magnetosphere, while changes in the dis- 
tribution parameters may be associated with changes in the way 
the magnetosphere- plasma interact i on pro ceeds. 

As discussed by iBurnard et al.l (Il983l) . plasma generally en- 
ters the magnetosphere of accreting pulsars according to vari- 
ous instabilities. These authors also conclude that, for the ob- 
served luminosities and spin-periods typical of bright accreting 
pulsars, the plasma mainly penetrates the magnetospheric bound- 
ary via Rayleigh-Taylor instabilities. If the accretion rate de- 
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Table 1. Best-fit parameters of the normal and "off-states" spectra with uncertainties at Icr confidence. 




I I I I I I I 



1.0 

Time. MJD-54634 




1.75 



Fig. 1. Observation-long lightcurve in the 0.4-1 2 keV energy range using data from all XIS units is shown in the bottom panel. The 
upper panels show close-up views of the three detected off-states. Here, the time axis is ticked every pulse period. 
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Fig. 2. Pulse profiles in four different energy ranges are shown for the normal (upper panels) and the "off-state" (bottom panels). 



creases, the rotating magnetosphere will inhibit accretion via 
Rayleigh-Taylor, therefore, for low-luminosity pulsars with in- 
termediate rotation r ates, the KHI accretion channel dominates 
(lBurnardetal.lll983h . 

The ways through which the plasma can penetr ate the mag- 
netosp here have been reviewed more recently by iBozzo et all 
(l2008h . who also provide estimates for the leak rates of various 
mechanisms. For a system with parameters similar to Vela X-1 
in "off-state", the hi ghest rate is e xpected to be provided by KHI 
(see section 3.2.2 of iBozzo et al.l for the details). The accretion 
luminosity is estimated in this case to be 

i'KH - GMnsMkh/^ns = 

7.4 X 10^^577KH<io(l + 16^Gio/(5/?Mio))'^'7^^^ erg s"' 

1 + Pi /Pe 

Here /?gio and /?mio are the capture and magnetosphere radius 
respectively, in units of 10'" cm; pi.e are the densities within and 
outside of the magnetosphere. According to Bozzo et al. (200^, 
77kh ~ 0.1 and the density ratio is estimated to be between 

i+Pi/Pe \0.1;7kh/J"'^mioV8 

where h is the fractional height of the area where the plasma 
and the magnetic field coexist, in units of the total thickness of 



the KHI unstable layer (iBurnardet al.lll983h . and -Ps283.5 is the 
spin-period in units of 283.5 s. We assume a canonical neutron 
star radius of /?ns ~ 10 km. In the case of Vela X-1, for the ob- 
served "off-state" luminosity of ~ 2.4 x lO-'^ erg s"', a magnetic 
field of B > 2 X 10'^^ G is required, if the KHI unstable layer is 
rel atively thin (h ~ 0.05 ), or B ~ lO'^'G, if h ~ 1, as suggested 
bv lBurnardetan(ll983l) . 

Evidence of such a high magnetic field in Vela X-1, sur- 
prisingly stronger than the one estimated from the CR SF energy, 
are extensively discussed in lDoroshenko et al.l(l2010bl) . Here, we 
wish only to point out that the discrepancy still can be cleared 
if the line formation region is located several kilometers above 
the neutron star surface. In fact, for the average observed lu- 
minosity of the normal state, ~ 4 x 10^^ ergs"', one expects 
that an accretion column with height up to ~ 10 km will arise 
(Lyubarskv & Sunyaev 1988; Doroshenko et al. 2010a). This im- 
plies a factor of ten decrease in field strength at top of the column 
given that B ~ -So((^ns + 10)//?ns)"'' ~ O.IBq. It is therefore suf- 
ficient to assume that the accretion column exists and that the ob- 
served CRSF forms closer to the top of the column to reconcile 
the strong magnetic field required by the observed "off-states" 
luminosity and the measured CRSF centroid energy. 

The observed change in the high-energy pulse profile is ex- 
plained under the assumption that an accretion column does in- 
deed exist. The sharp dip, evident at high luminosity and softer 
energies, would then be due to the eclipse of the polar cap by 
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Fig. 3. Histogram of the measured count-rate of the XIS 
lightcurve (solid line). The same is shown for the off-states 
(shaded) and for the rest of the lightcurve separately. 

the accretion column. As the luminosity drops and the column 
ceases to exist, the hard X-rays can pierce through and the po- 
lar cap is observed directly: a pronounced peak is observed in- 
stead of a dip. In other words, the high-amplitude peak, which 
appears in the "off-state" pulse profile at hard energies around 
pulse phase 0.75 can be attributed to the direct emission from 
the polar cap. In this scenario, the accretion stream would still 
absorb the soft X-rays so the dip is still observed at lower ener - 
gies. A similar scenario is discussed bv iKlochkov et al.l (l2008l) 
to explain the pulse profile variations of EXO 2030-1-375 during 
outbursts. 

4. Conclusions 

We presented an analysis of the "off-states" of Vela X-1 ob- 
served during a 100 ks Suzaku observation of the source. For the 
first time ever, we have been able to study and characterize the 
properties of Vela X-1 during the so-called "off-states", detect- 
ing pulsations and measuring the spectrum and flux distribution. 
We also confirm the presence of a CRSF at 25 keV in the normal 
state of the source. We conclude that these observational results 
strongly suggest that the emission is still powered by accretion 
and that the drop in luminosity has a magnetospheric origin. The 
observed X-ray luminosity of the "off-state" may be naturally 
explained in the gate d accretion scenario originally proposed by 
Bumard et al.l (Il983l) and recently investigated by (iBozzo et al.l 
2008h . if the neutron star is strongly magnetized and the plasma 
enters the magnetosphere via Kelvin-Helmholz instabilities. To 
our knowledge this is the first time that the theoretical consid- 
erations of magnetosphere-plasma intera ctions proposed so far 
(iBurnard et al.l 119831: IBozzo et al.1 l2008h find an observational 
confirmation. Same interpretation might apply to several other 
sources for which episodes simil ar to the "off'-states" observed 
in Ve la X-1 have been reported (iNaik et al.ll2005t iGogii^ et al.l 
l20Tlh . 
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